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Abstract 

We investigate a class of localized, stationary, particular numerical solu- 
tions to the Maxwell-Dirac system of classical nonlinear field equations. The 
solutions are discrete energy eigenstates bound predominantly by the self- 
produced electric field. 
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I. INTRODUCTION 



There are many examples of classical solitary wave solutions to nonlinear field theory 
equations. Some of these are useful in quantum field theory as a stationary point in the 
action functional about which one quantizes the field Although these solutions are 
often relegated to model equations in fewer than three space dimensions 0, we consider the 
Maxwell-Dirac system of equations in 3+1 dimensional space-time, a nonlinear system of 
PDE's involving twelve real functions of four variables. The general solution of this system 
is certainly well beyond our grasp; however, we will obtain a class of particular solutions by 
making simplifying assumptions and utilizing numerical methods. 

Before we embark on a search for a localized solution it would be wise to consider on 
what grounds such a solution is plausible. It is well known that when one considers the 
Dirac equation with an external "repulsive" potential the possibility arises to obtain bound 
state solutions 0. This potential produces bound states with discrete energies that rise from 
the continuum of free negative energy states in the same way that an "attractive" potential 
produces lower energy states from those of positive energy. For a great enough "repulsive" 
potential we may even obtain bound states of positive energy. This interesting phenomenon 
goes under the name of "Klein's paradox" and provides our motivation. In our case the 
"repulsive" potential is provided by the charge feeling its own electric field. 



II. FIELD EQUATIONS 

The Maxwell-Dirac equations are obtained from the Lagrangian density 

L = i^-i^'d^^ - fr^ _ q^l^^A^ - ^F^'F^,, (2.1) 

in which the c-number fields are the Dirac spinor \Ef, which can be considered a four com- 
ponent single-column matrix, and we have = A), # = \E'"I'7°, F^'^ = d^A'^ — d^A^, 
dfj_ = = (t, x), and 7^ are the 4x4 matrices 
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7° = I 1,7 



7 \ / <j\ 

.0 -// V-^ / ' 

in which / is the 2x2 identity matrix and o represents the Pauh matrices 

1\ /O -i\ l\ 

1 o; ; Vo -1 

The Euler-Lagrange equations apphed to Eq. ( p.l|) give the Maxwell-Dirac equations 

7^(29^ - qA^)^ - ^ = 0, (2.2) 

d^F^"" = q^^f"^. (2.3) 

Throughout we use natural units, in which we have rescaled length, mass, and time so that 
h = m = c = 1. 

We now begin making assumptions about the solution we wish to look for. We require 
that satisfy the Lorentz condition d^A^ = and that \E' is an energy eigenstate, = 



i){x)e-'^K Equation now reduces to Poisson's equation 

V^A'^ = (2.4) 

where is the 4-current 

f = q^Y^^ (2.5) 

and Eq. ( p.2|) can be written as a Hamiltonian eigenvalue equation 

Etp = Htp= [7°7 • (-iV - qA) + 7° + q^ip. (2.6) 

We would now like to assume spherical symmetry for our wave function; however, we 
find that the resulting vector potential A has an angular dependence that destroys the 
symmetry. This gives us two options: we may throw out the magnetic term in the hope that 
its contribution will be small and solve the resulting one dimensional, spherically symmetric 
problem, or we may attack the non-spherical problem. We proceed with the former and 



save the latter for Sec. [V|. By removing A we reduce the Maxwell-Dirac system ( [2.1D to a 
massless scalar-Dirac system 



(2.7) 



Note that, with 7" in the couphng term, our massless scalar-Dirac system ( p.7[ ) is not Lorentz 
invariant and is therefore only of use in approximating the Maxwell-Dirac system. 

The necessity for a spherical charge distribution restricts our wavefunction to four pos- 
sible configurations corresponding to total angular momentum up or down, = ±|, and 
the quantum number k = ±1, 
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—g{r) cos9 


in spherical coordinates {r,9,(j)), where we have chosen to orient the angular momentum 
along the z-axis. Note that the two k = —1 configurations have the same angular dependence 
as the Hydrogen ground state. Equations ( |2.4| - |2.(j| ), neglecting A, now reduce to the radial 
equations 



$" + -$' = -g(/2 + /). (2.10) 
r 



We symmetrize the Hamiltonian in Eq. ( p. 91) by the similarity transformation F = r/, G = 
r^f, and further simphfy Eqs. ( p.9|) and ( p.lOj ) by identifying the fine structure constant 
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a = 1^ and defining the potential V = q^, giving us 



^rV)" = -^{F^ + G^). (2.12) 
r 

We also restrict our wavefunction by imposing the normalization condition 

//■oo 
^3^^t^ = A'K j dr{F^ + G^), (2.13) 



ensuring a total charge of q. Note that niz does not appear in Eqs. ( |2.11[ - PTTBD ; hence the 
energy levels are independent of the choice m^. = 



III. SOLUTION 



We may readily solve our system of O.D.E.'s ( |2.11| , |2.12|) by a variety of means, including 



power series solution 0], Pade series approximation ||^, and numerical methods. We choose 
the latter as it seems the most straightforward path to obtaining a solution. 

We discretize Eqs. ( p.ll| ) and (|2.12|) using second order finite differences to obtain a 



standard, linear, symmetric, matrix eigenvalue problem for {G{r), F{r)) and E coupled 
nonlinearly to a symmetric, matrix inverse problem for rV. We obtain solutions by solving 
each linear problem independently and iterating to convergence from an initial guess. This 
technique is similar to Newton's method of solving the nonlinear system but allows for better 
behavior in solving the eigenvalue problem at the cost of slower convergence. We use inverse 
iteration for the eigenvalue problem together with the conjugate gradient method for 
the matrix inverses. Although more efficient methods are available for the one dimensional 



problem [0, we have chosen methods which will work equally well for the case of two 
independent variables discussed in Sec. |V|. 

Figure |I] shows a solution to our system ( |2.11[ - |2TT3D for the choice of a = 2.7. We may 
use the numerical solution to calculate expectation values, such as (r) and the strength of 
the neglected magnetic interaction. Table | shows the results for various selections of the 
parameter a, which completely determines the set of solutions. The energy E as a function 
of a, K, and the number of nodes in g, represented by n, is shown in Fig. ^. We can see from 
Fig. H and directly from Eq. (|2.11| ) that localized solutions are only possible for —1 < E < 1. 
Note that positive energy solutions exist for a > 2.4 and that there is a spectrum of large n 
states of energies near negative one for any choice of a. 



IV. MAGNETIC INTERACTION PERTURBATION 

For the solution to our scalar-Dirac equations to be a viable approximate solution to the 
Maxwell-Dirac equations, we must establish that the inclusion of the magnetic interaction 
does not significantly alter the solution. We do this properly in Sec. by modeling the larger 
system ( |2.4| - |2.6| ), including the angular dependence; however, we will use perturbation theory 
with our spherically symmetric solution to determine the approximate energy shift AEm 
due to the magnetic interaction and find good agreement with the full Maxwell-Dirac case. 

The magnetic interaction term in the Hamiltonian density is 

Hm = -3-A. (4.1) 

We use our solution form ( ^.8| ) in Eq. ( [^.5| ) to get the current 

j = 4qKmzf{r)g{r) sin 9(j). (4.2) 

From this we obtain the vector potential A = A{r) sin Ocj), in which A is calculated via the 
Green's function integral 

A{r) = - / dnrl^jin), (4.3) 

o Jo Ty 
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where r< (r>) represents the lesser (greater) of r and ri. We may now calculate the approx- 
imate magnetic energy shift to first order via 

AEm ^ J d\HM (4.4) 
128 o /""^ /""^ , r< 



i fOO POO f 

-Ti^a / drF{r)G{r) / dr^^F{ri)G{ri). 
Jo Jo r^^ 
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The values obtained for AEm are seen in Fig. ^ and Table H and found to be small compared 
to the binding energy. 

We may also use our current ( [4.2|) to calculate the magnetic moment via 

f^=\j X /) (4.5) 

= 2qKmziT'^ / drrF{r)G{r)z. 



Note in Table | that we obtain the result — f when we choose n = 1. 



V. FULL MAXWELL-DIRAC SOLUTION 



From our calculations in the previous section we expect that the inclusion of the magnetic 
interaction will result in the emergence of a small non-trivial angular dependence on the 6 
coordinate. We also expect the solution to maintain its axial symmetry and proceed in 
cylindrical coordinates (p, z, (p), assuming a wavefunction of the form 

iIj2{p, 2)e^(™-+5)<^ 

\-i^i{p, 2)e^(™^+5)'^ ) 

and potentials of the form $ = ^V{p, z),A= ^A{p, z)cj). Equations ( |2.4| - |275| ) now reduce to 



(5.1) 



1 



and the Hamiltonian (12. 61) becomes 



(5.2) 
(5.3) 
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Note that our wavefunctions in both cases ( |2.8| ) and ( ^.1| ) are eigenfunctions of the z com- 
ponent of the angular momentum with eigenvalue = ±i, where 



J z = Lz -\- l^zi 



in which L = —i{x x V) and 




(5.6) 



(5.7) 



The full system ( ^.2[ - ^.5|) is once again spin degenerate since the change to = — m'^ 
produces the same system with the change {A = — A', \E'i = — ^^2' ^2 = ^'n ^3 = ^4 = 
— \E'3). Although our scalar-Dirac wavefunctions (|2.8| ) are eigenfunctions of K = '-/^{a- L + 1) 
with eigenvalue k, our solutions of the non-spherically symmetric case will vary slightly from 
these eigenfunctions. 

Figure |^ shows a numerical solution to the full system for the same parameter choice as 
in Fig. Note that the angular dependence is virtually indistinguishable from the approx- 
imate solution and that the energy agrees reasonably well with the first order perturbation 
approximation. 



VI. CONCLUSION 



We have found a class of solutions to the full Maxwell-Dirac equations and good approx- 
imate solutions via a scalar-Dirac equation. In practice, the approximate solution is easier 
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to work with and provides better accuracy for most calculations. 

The issue of stability has not been directly addressed. The success of our iterative 
solution method suggests that each solution is stable in regard to slow collapse or expansion; 
however, we suspect that each solution will be unstable via radiative transitions to states of 
large negative energy, as is any bound state solution to the Dirac equation. 

The interpretation of our solution is not immediately clear. Several authors have used 
similar solutions to construct hadrons from interacting quarks ; however, these solutions 
were to systems of nonlinear scalar fields interacting with the Dirac field and the solutions 
were essentially a result of the nonlinear scalar self-interaction rather than the coupling term. 
As it exists now, we see that our class of solutions to the Maxwell-Dirac system may not 
be immediately interpreted as representing the leptons. This is clear from the large value 
of < r > and negative value for E for well as from the experimental fact that 

weak interactions play the starring role in lepton transitions. It is not inconceivable that 
our solution could represent the leptons if we rework our analysis to accommodate the weak 
interaction. Such an undertaking would present an interesting avenue for future work. 

ACKNOWLEDGMENTS 

1 would like to express my appreciation to Roger Dashen and Henry Abarbanel for their 
encouragement and guidance and to the Hopgood Foundation and the San Diego chapter of 
the ARCS Foundation for their support. 



9 



REFERENCES 

[1] R. Rajaraman, Solitons and Instantons (North-Holland, Amsterdam, 1982). 
[2] A. Das, J. Math. Phys. 34, 3986 (1993). 

[3] J. J. Sakurai, Advanced Quantum Mechanics (Addison- Wesley, N.Y., 1967). 
[4] M. E. Rose, Relativistic Electron Theory (John Wiley and Sons, N.Y., 1961). 
[5] W. L. Ditto and T. J. Pickett, J. Math. Phys. 29, 1761 (1988). 
[6] W. Press, Numerical Recipes in C (Cambridge University Press, N.Y., 1992). 
[7] S. Salomonson and P. Oster, Phys. Rev. A 40, 5548 (1989). 

[8] W. L. Bardeen, M. S. Drell, M. Weinstein, and T. M. Yan, Phys. Rev. D 11, 1094 (1975). 
[9] T. D. Lee and G. C. Wick, Phys. Rev. D 9, 2291 (1974). 



10 



FIGURES 

FIG. 1. A normalized, localized solution to our scalar-Dirac system for the choice of a = 2.7, 
K = 1, and n = 0, showing /(r), g{r), and V{r) scaled down by a factor of 6 to be fully visible. 
Note the asymptotic behavior of V{r) ~ ^ for large r. The expectation value < r > is also shown. 
All values are in natural units. The one dimensional mesh was discretized into 200 points for this 
calculation. 

FIG. 2. Energy levels as a function of a for our scalar-Dirac approximation. The dashed curves 
represent deviations due to the magnetic interaction calculated via the first order perturbation. 
All values shown are in natural units. 

FIG. 3. A normalized, localized solution to the Maxwell-Dirac system for the choice of a = 2.7, 
showing the four components of the wavefunction V'i(P) z), the potential V{p, z), and the (f) compo- 
nent of the vector potential A{p, z). The contours go from min(light) to max(dark) and all values 
are in natural units. The two dimensional mesh was discretized into 60 x (2 x 60 + 1) points for 
this calculation. The energy was calculated as E2 = 0.22, to be compared with that obtained from 
the 1-D approximation, E + AEm = 0.30. 
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TABLES 

TABLE L Numerical results for our scalar-Dirac approximation on a 200 point mesh for several 
choices of a, n, and k. All values shown are in natural units. 



a 


n 


K 


E 


AEm 


< r > 


f^z/iqrri;,) 


1/137 





1 


-1 + 9 X 10-*^ 


-0.00 


572 


-1.2 


0.1 





1 


-1 + 2 X 10-^ 


-0.00 


42 


-1.2 


1.0 





1 


-0.83 


-0.00 


4.1 


-1.1 


2.0 





1 


-0.28 


-0.05 


1.9 


-1.0 


2.4 





1 


0.09 


-0.10 


1.5 


-0.9 


2.7^ 





1 


0.46 


-0.16 


1.2 


-0.8 


3.0 





1 


0.93 


-0.24 


1.0 


-0.7 


1/137 


1 


-1 


-1 + 3 X 10-6 


-0.00 


1710 


-0.4 


2.0 


1 


-1 


-0.75 


-0.00 


5.4 


-0.3 


3.4 


1 


-1 


0.06 


-0.06 


2.0 


-0.2 


3.8 


1 


-1 


0.89 


-0.14 


1.2 


-0.1 


3.0 


1 


1 


-0.68 


-0.00 


6.2 


-1.1 


5.0 


1 


1 


0.32 


-0.05 


2.4 


-0.9 


5.5 


1 


1 


0.90 


-0.10 


1.8 


-0.8 


2.0 


2 


-1 


-0.93 


-0.00 


18 


-0.4 


5.8 


2 


-1 


0.02 


-0.03 


3.5 


-0.2 


6.0 


2 


-1 


0.18 


-0.03 


3.1 


-0.2 


6.0 


2 


1 


-0.37 


-0.01 


6.2 


-1.0 


6.0 


3 


1 


-0.71 


-0.00 


13 


-1.1 



^Solution shown in Fig. ||. 
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This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-th/9410244v2 
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